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Abstract 

Supermassive black holes (SMBH) are typically surrounded by a dense stellar population in galactic 
nuclei. Stars crossing the line of site in active galactic nuclei (AGN) produce a characteristic transit 
lightcurve, just like extrasolar planets do when they transit their host star. We examine the possibility 
of finding such AGN transits in deep optical, UV, and X-ray surveys. We calculate transit lightcurves 
using the Novikov-Thorne thin accretion disk model, including general relatistic effects. Based on the 
expected properties of stellar cusps, we find that around 10® solar mass SMBHs, transits of red giants 
are most common for stars on close orbits with transit durations of a few weeks and orbital periods 
of a few years. We find that detecting AGN transits requires repeated observations of thousands of 
low mass AGNs to 1% photometric accuracy in optical, or ~ 10% in UV bands or soft X-ray. It may 
be possible to identify stellar transits in the Pan-STARRS and LSST optical and the eROSITA X-ray 
surveys. Such observations could be used to constrain black hole mass, spin, inclination and accretion 
rate. Transit rates and durations could give valuable information on the circumnuclear stellar clusters 
as well. Transit lightcurves could be used to image accretion disks with unprecedented resolution, 
allowing to resolve the SMBH silhouette in distant AGNs. 

Subject headings: accretion disks — galaxies: active — quasars — techniques: photometric 



1. INTRODUCTION 

Photometric transits have been successfully used to 
identify and characterize transiting extrasolar plan- 
ets since the discovery o f the first one, HP 209458b 
(jCharbonneau et all 120001: iHenrv et al.1 120000 . Transit 
shape is a telltale of planetary, orbital, and stellar pa- 
rameters. Moreover, the apparent time-dependent red- 
shift of the star due to the planet covering its receding or 
approaching side during transit can reveal the projected 
angle between the planetary orbital axis and the stel- 
lar rotational axis (Ros siter-McLaughlin effect, Rossiterj 
119241: iMcLaugh lin 19241. These methods show that tran- 
sits are very powerful in probing planetary systems. 

Similarly, active galactic nuclei (AGN) accretion disks 
can be probed by observing occultations in X-ray due to 
broad line region clouds: optically thick clouds orbiting 
the supermassive black hole (SMBH) in the region from 
where broad emission lines of the AGN are thought to 
originate. These occultatio ns have a large covering fac- 
tor of ~ 0.1 to 1 (see e.g. iMcKerna n fc Yaqoo bI Il998l: 
Turiier et aP [2 008: Bianchi"!eta]J |2009; Maiolin o et al 
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( 2011aD pointed out that analogously to the Rossiter 



McLaughlin effect, temporally resolved spectroscopic ob- 
servations of an eclipse could be used to probe the appar- 
ent temperature structure of the accretion disk and the 
origin of the iron Ka emission line, allowing to constrain 
the black hole spin and inclination. 

In this paper, we examine the possibility of stars on 
close orbits transiting their host AGN. There are mul- 
tiple coincidences that make it possible to detect such 
events: the radius of large stars is comparable to the 
characteristic size of an accretion disk around a 10® Mq 
SMBH, resulting in transits deep enough to detect. The 
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orbital period is a few years in the innermost regions 
where the stellar number density is highest, short enough 
for repeated observations. Finally, the typical transit 
duration for these innermost stars is an hour to a few 
weeks (depending on the stellar population and the ob- 
serving band), which is longer than the typical cadence 
required to observe AGNs, but still accessible on human 
timescales. 

In order to produce a detectable signature in the 
lightcurve, the transiting object has to be a main se- 
quence O or B star, a Wolf-Rayet (WR) star, an AGB 
star with a surrounding dust cloud (OH/IR star), a 
young supermassive star, or an evolved giant. Late type 
main sequence stars are too small to cause a transit de- 
tectable from the ground in an AGN with black hole mass 
> 10^ Mp, , unless they are "bloated" by irradiation of the 
AGN (|Lohfink et al.ll20Tl . 

High photometric and temporal resolution observa- 
tions of AGN stellar transits have the potential to map 
the accretion disk structure with an unprecedented ac- 
curacy. Stars are optically thick in all electromag- 
netic bands, and un like broad line region clouds (e.g. 
iMaiolino et al.l[2010( ) , they have a simple spherical geom- 
etry, making it easier to interpret lightcurves and to re- 
construct the image of the accretion disk along the transit 
chord (projected stellar trajectory). Furthermore, such 
transits offer unique observations of individual stars in 
distant galaxies. Detections of multiple transits could re- 
veal valuable information on the number density of stellar 
cusp central regions. 

In this paper, we calculate transit depths, durations, 
periods, rates, and instantaneous probabilities based on 
physical models of nuclear stellar clusters. We also 
present simulated transit lightcurves and transit depth 
maps in multiple frequency bands. The shape of the 
lightcurve depends on the observing wavelength, the 
mass and spin of the SMBH, the accretion rate, the inch- 
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nation of the accretion disk, the projected position of the 
transit chord, and the orbital velocity of the transiting 
object. Therefore observations with sufficient photomet- 
ric accuracy and time resolution allow us to constrain 
these parameters, and to test the employed accretion and 
general relativity models. 

In Section [21 we review observations of large stars 
closely orbiting Sgr A*, the SMBH at the center of the 
Galaxy ( ij2.ip . and summarize theoretical models of semi- 
major axis distribution and mass segregation ( §2.2p . We 
set up models of the stellar population and radial dis- 
tribution ( ij2.3p . and determine the minimum (i j2.4[) and 
maximum fi j2.5l) orbital radii for transits. We present 
accretion disk thermal emission models and ray-tracing 
simulations in Section |3l showing transit spectra fi )3.ip . 
transit depth maps fi J3.2l) and lightcurves (t j3.3|) . We cal- 
culate transit probabilities in Sectional In Section [SJ we 
determine the local density of AGNs of interest fi jS.ip . 
and study the feasibility of transit detections with opti- 
cal m5.2[ H5.3^ and X-ray instruments f ii5.4l) . The most 
important simplifying assumptions, caveats, and impli- 
cations are discussed in Section [B] Finally, we conclude 
our findings in Section [T] 

2. NUCLEAR STELLAR CLUSTERS 

In this section we review observations and models of 
stellar distribution in galactic nuclei. 

Let MsMBH denote the mass of the SMBH, and define 
Mg = Msmbh/(10^ -^0) and the gravitational radius 
Rg = GMsmbh/c^. Then 

i?g 4.8 X 10"^ mpc Me = 2.1 Rq Mq , (1) 

where 1 mpc — 10^'^ pc — 206 AU. For non-spinning 
black holes, the Schwarzschild radius is i?s = 2i?g, 
and the radius of the innermost stable circular orbit 
is -Risco — GRg- For maximally spinning black holes, 
Rs — Rg, and i?isco — Rg for prograde orbits in the 
equatorial plane. 

2.1. Stellar populations 

Many ga l axies h ost a dense stellar cusp in their nucleus. 
ISeth et al.l ()2008aD find that galaxies with a massive nu- 
clear cluster are more likely to be active. Stars captured 
and transported inwards by the accretion disk may serve 
to fuel the AGN (Miralda-Escudc & Kollmcicr 20Q1). 

In the Galaxy, lEisenhauer et al.l (|2005i) report the or- 
bital parameters of six B type main sequence stars or- 
biting the central SMBH on highly eccentric orbits with 
semi-major axes less than 16 mpc ~ 10^ Rg. These stars 
may be remnants of stellar binari es tid a lly dis rupted by 
the SMBH, as first p roposed by iHillsl (| 19881) (see also 
lYu fc Tremainel[2003D . Candidates for the binary coun- 
ter parts ejected with high velocity have been identified 
b ylBrown et al.l (I2009D . 

iBartko et al.l ()2009[ ) find more than a hundred O and 
WR stars further from the Galactic Center, within 
1 pc = 5 X 10^ Rg projected distance. These stars could 
have formed in a massive s elf-gravitating gaseous disk 
(e.g. lLe^[2007L iHobbs fc Navak shin 2009), or fo rmed 
furthe r away and captured in close orbits by the IHillsl 
(|1988D mechanism or by a clus ter of stellar mass black 
holes ([Alexander fc Livioll2004[ ). Young stars could also 
be delivered to this region by an infalling globular cluster 
(fTremaine et al...l975; .Milosavlicvic fc Mcrritt.,2001i) . 



Most main sequence stars in the vicinity of a SMBH 
eventually evolve into red giants or supergiants, our can- 
didates for transiting the AGN. Many such giants have 
been observed within 1 pc of the Galact ic Center (e.g. 
'Genzel et ^ 120101: IPaumard et"all 120061: IBartko et all 
i2009h .20101 ). The fraction of stars in the giant phase 
within a stellar population depends strongly on its ini- 
tial mass function (IMF) and formation history. 

A more exotic possibility is the formation of ~ 10^ Mq 
supermassive stars in the accretion disk as suggested by 
[Goodman fc Tan (2004). Such a star would form outside 
1000 Rs, have a radius of approximately 360 Rq in case 
of solar metallicity, and would radiate at its Eddington 
limit, being luminous enough to have a detectable optical 
photometric signature not only when it transits the AGN 
but also when it is occulted behind it. 

Observations of the Galactic nucleus show that the 
innermost cluster of young stars (so-called S- stars) 
is isotropically distributed ([Genzel et al.l l2010f ). as 
predicted by theoretical arguments. Even if stars 
are forme d on an orbit co p lanar with the accre- 
tion disk. Ra nch &: Tremaini ([1996*) provide a relax- 
ation mechanism that could rapidly randomize the or- 
bital orientations. The possible presence of interme- 
diate mass black holes may help catalyze this process 
([Gualandris fc Merrittl [2009) . And even if stars cluster 
in disks, this coherent behaviour averages out when ob- 
serving multiple galactic cores as long as the stellar disks 
and the accretion disks have independent orientations. 
Therefore we assume an isotropic distribution of stellar 
orbits in the cluster for the purpose of our probability 
estimates. 

2.2. Density profile 

We assume circular orbits for simplicity, and denote 
orbital radius by r. A star on an eccentric orbit with 
the same semi-major axis would produce a transit of the 
same depth, with a transit probability and transit length 
w ithin a factor of order unity. 

iBahcall fc WoU (!l976) showed that the equilibrium 
spatial number density distribution of a stellar cluster 
around SMBH is proportional to r~" with a = 1.75 if 
the stars in the cluster have the same mass. Analyti- 
cal and numerical investigations of multiple mass pop- 
ulations show that the distribution for each mass bin 
is still likely to follow a power law. The value of a 
is predicted to be smaller (down to 1.3) for lower 
mass stars ([Ba hcall fc Wolj Il977t iFreitag et aP 120061 : 
iHopman fc Alexander. .2006) . For massive stars repre- 
senting a small mass fra ction in the stellar cluster, a can 
be between 2 and 2 .75 f Alexand er fc Hopmanll2009D . or 
as large as 3 (IKeshe t et al. 2009). Observations of the 
Galaxy by IBartko et al.l ([2010) show that WR/0 stars 
from a distance of 30 mpc out to 600 mpc form an inter- 
esting anisotropic structure called the clockwise disk in 
the Galacit nucleus. These stars are distributed with a 
density exponent a — 2.4 ± 0.2, while the B star popu- 
lation from 30 mpc to 1 pc exhibits a — 2.5 ± 0.2. Note, 
however, that the age of main sequence O stars and WR 
stars is less than the tw o-body relaxation ti me ~ 0.1- 
1 Gyr at r - 100 mpc ([O'Learv et al.ll2009( ). therefore 
they are not expected to represent the equilibrium state. 
The observed mass distribution of solar mass stars in 
the Galactic nucleus is fit by a broken power law with 
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a = 1. 2 and 1.75 inside a nd outside of 0.22 pc, respec- 
tively (|Sch6del et all 120071) . 

To estimate the total number of stars, we first define 
the radius of influence r; as the radius of the sphere cen- 
tered on the SMBH containing a total mass of 2 A/g]y[g]^ in 
stars and stellar remnants. In case of a singular isother- 
mal sphere (a = 2), this equals to GMsmbh/ct^, where a 
is the velocity dispersion of stars further than ri (jMerrittl 
|2004[) . In order to get an estimate of the number of stars, 
we set n — GMsmbh/c^, independently of a. 

Using the M-a relation now allows us to express the 
radius of influence as a function the supermassive black 
hole mass only, in the form 



n = ro M7 . 



(2) 



Here ro and 7 depend on the coefficients of the M-a 
relation. In particular, 7 = 1 — 2 //?, where is the slop e 
of logM-logcr, as defined b y e.g.pTremaine et al.l (l2002| ). 
For example, the best fit of iTremaine et al.l (|2002| ) (with 
P = 4.02) results in ro = 1.234 pc an d 7 = 0.50, and 
the best fit of 'Fcrrarcsc & Ford (I2005D (with /3 = 4.86) 
results in rp — 0.881 pc and 7 = 0.59. Fo r our numerical 
results, we adopt the best fit values of iGiiltekin et al.l 
(|2009) (with 13 = 4.24): rp = 1.075 pc and 7 = 0.53 (but 
use 1 pc to normalize rp in our parametric expressions). 

Now let us consider the stellar population in the vicin- 
ity of the SMBH. We assume that this population con- 
tains a species of stars with mass A/^,, having a radius 
Ri, large enough to produce detectable transits. We also 
assume that the number density of these stars follows the 
power law r^". Let b denote the mass fraction of these 
large stars within the sphere of influence relative to the 
total mass of all stars and stellar remnants. 

Typically the inner cutoff radius rmin for the stellar dis- 
tribution is much smaller than the radius of influence (see 
i) l2.4l for numerical justification). Under this assumption, 
the spatial number density of the stars in question as a 
function of orbital radius is 



n{r) — b 



3 - a Af; 



SMBH r 



2n A'h 



(3) 



Note that this argument has two weaknesses: first, the 
M-a relation has a relatively large scatter. For a given 
SMBH mass, the intrins ic scatter of the bulg e velocity 
dispersion is ~ 0.3 dex (jGiiltekin et al.l [20091 ). Second, 
we applied results for the isothermal sphere to power 
law distributions with different exponents. This limits 
the accuracy of the transit rate estimates presented in 
Section m 

2.3. Adopted models 

We consider two simple models for the transiting stellar 
populations around AGNs, summarized in Table [TJ 

First, we assume a young stellar population, motivated 
by the observations of a young population with an ex- 
tremely top-heavy initial mass function of mean stellar 
mass ~ 30 A/q following a density pro file with a » 2. 5 
in the central parsec of the Galaxy (jBartko et al.lf20Tol ). 
Thus in our model, we assume that a fraction of the 
stars are O type, with stellar mass Mi, = 30 Mq, radius 
i?^, = 11 i?0, and a = 2.5. This exponent is consistent 
with observations and theoretical predictions reviewed in 
Section [2?2] We assume that these O stars constitute a 



mass fraction h = 0.01 of the population. This is consis- 
tent with the estimated total mass of WR/0 stars in the 
Galactic Center if we consider that these stars are con- 
fined in the center part of the sphere of influence. How- 
ever, note that this mass fraction depends sensitively on 
recent star formation rate and initial mass function of 
stars in the neighborhood of the SMBH. 

For the second model, we consider an evolved, relaxed 
population of stars, and optimistically assume that a 
b — 0.01 mass fraction of them are giants (or main se- 
quence stars otherwise enlarged, like "bloated" or sur- 
rounded by a dust cloud), which are large enough to 
produce detectable AGN transits. For these giants, we 
assume i?^, = 110 i?0, A^* = 1.5 Mq, and a Bahcall-Wolf 
equilibrium value of a = 1.75 (see Section [2?2|) . 

Note that depending on the star formation history, O 
stars and red giants might coexist in the cusp, in which 
case their contributions to transits add up. 

2.4. Minimum orbital radius 

The inner orbital radius cutoff of the stellar distribu- 
tion, Tinin, is set by gravitational wave inspiral, and tidal 
and coUisional disruption. While gravitational w ave in- 
spira l is the limiting factor for compact objects (jPetersI 
Il964| ). tidal or coUisional disruption sets a tighter con- 
straint for stars. 

The tidal disruption radius is 



2 -^^SMBH \ 



0.013 mpc Mi 



Af. 
30 M, 



o 



R. 



11 Rr. 



(4) 



where rj ranges from 0.8 to 3.1 depending on the com- 
press ibility and polytrop ic mass distribution index of the 
star (jPiener et al.lll995f ). We adopt the moderate value 
77 = 2. The tidal disruption radius is given here normal- 
ized for a massive main sequence star. 

CoUisional disruption might be responsible for the de- 
pletion of luminous gian t stars in the i nner 80 mpc 
of the Galacti c Cen ter ([Alexanderi I1999D . Following 
IHopman et al.l H2007f ). one can write that the rate at 
which collisions decrease stellar density is 



dn{r,t) n'^{r,t)v{r)Y. 



dt 



A^coll 



(5) 



where E = irR^ is the coUisional cross-section, and on av- 
erage, A^coU ~ 30 co llisions are required to disrupt a star 
(jFreitag et al.|[2006l ). We define the radius limit for coUi- 
sional disruption as the radius where the stellar density 
e-folds in time t if the initial value is given by Equation 

m 



''coll = ro 



1.33 X 10"^ 6(3-a)Afj 



10^ yr / \ 1 pc 



ro 



AT. 
30 Mp 



i-7(3-Q) 



R. 

11 Rr. 



(6) 



We set the timescale to be t = 10^ yr: this is within an 
order of magnitude of both the main sequence lifetime 
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Table 1 

Orbital radius and period limits for AGN transits for the adopted stellar population models (b = 0.01 for each model) 



MSMBH 


model 








a 




»'coll 


^min 


^Iciis 






^orb C'^min) 


^orb (^max) 


(Mq) 




(Re) 




(Mq) 




(mpc) 


(mpc) 


(mpc) 


(mpc) 


(mpc) 


(mpc) 


(y) 


(yr) 


10^ 


O stars 


11 


51.8 


30 


2.5 


0.006 


0.16 


0.16 


11 


320 


11 


0.57 


340 


10^ 


red giants 


110 


518 


1.5 


1.75 


0.16 


0.54 


0.54 


22 000 


320 


320 


3.7 


53 000 


10*^ 


O stars 


11 


5.18 


30 


2.5 


0.013 


0.40 


0.40 


11 


1075 


11 


0.75 


110 


10*^ 


red giants 


110 


51.8 


1.5 


1.75 


0.35 


1.3 


1.3 


22 000 


1075 


1075 


4.2 


100 000 


10^ 


O stars 


11 


0.518 


30 


2.5 


0.028 


1.0 


1.0 


11 


3 600 


11 


1.0 


34 


107 


red giants 


110 


5.18 


1.5 


1.75 


0.76 


3.0 


3.0 


22 000 


3 600 


3 600 


4.9 


210 000 



of early type stars, and the lifetime of the giant phase 
for less massive stars. Note that here we only consider 
collisions within the large species, not with other, much 
smaller stars, which are less likely. 
The minimum orbital radius is 

fmin = niax(rt,rcoii) . (7) 

Table [U lists the minimum and maximum orbital radii for 
the two stellar species and three different SMBH masses. 
We find that in every case, collisions set a tighter con- 
straint than tidal disruption for the potentially transiting 
stars. 

Table [1] also lists the Keplerian orbital periods for stars 
at the minimum and maximum radii, which can be cal- 
culated as 

3 

To.. = Syr[^y M-K (8) 

We find that the closest main sequence stars have or- 
bital periods of approximately one year, making repeated 
transit observations feasible. However, collisions set a 
larger radius limit for giants, resulting in longer orbits. 
This inner radius limit depends on the number fraction 
of giants: smaller b implies less frequent collisions and 
thus allows closer orbits. However, a smaller value for b 
would also mean smaller transit probabilities for a given 
AGN (see Section H below). 

2.5. Maximum orbital radius 

For the maximum orbital radius r,„ax of stars capable 
of producing a transit, we have to consider two factors: 
gravitational microlensing due to the star, and the valid- 
ity range of the presumed number density power law. 

Gravitational microlensing caused by the transiting 
star can bend the light rays of the AG N which niay sig - 
nificantly distort the transit lightcurve ()Pacz vnskill 1 9860 . 
This can happen if the transiting object is farther from 
the AGN than the radius ricns at which the Einstein ra- 
dius equals the apparent angular radius of the transiting 
object: 

= '''''' iirk) (^) • 

Therefore we restrict our transit probability calculations 
to the contribution of stars within orbital radius ricns- 

Note that this is a different configuration than a galaxy 
microlensing a distant quasar, which can also be used 
to probe the spatial structure of accretion di sks around 
AGNs (e.g. for the case of Q2237+0305, see iKochanekl 
12004 and references therein). 



The power law distribution discussed above only ap- 
plies to the stellar populations within the radius of influ- 
ence from the SMBH, where its gravitational field dom- 
inates. In this paper, we do not investigate the distribu- 
tion of stars outside the sphere of infiuence, but conser- 
vatively ignore their contribution to transit probabilities. 

The maximum orbital radius is thus the smaller of the 
microlensing radius and the radius of infiuence: 

''max = min(nons, n) . (10) 

Table [1] shows that typically microlensing is the limiting 
factor among main sequence O stars, whereas the radius 
of infiuence limits transits of the much less dense red 
giants. 

3. TRANSIT OBSERVABLES: SPECTRA AND 
LIGHTCURVES 

Next we derive the AGN spectra and the transit ob- 
servables. 

The AGN luminosity is bound ed by the Eddington 
limit I/Edd = 1.3x1044 erg s'^ Mq (jShapiro fc Teukolskvl 
Il983f l. We assume an Eddington ra tio of . 25 as 
our fiducial value , following (Kollmci er et al.l 120061 : 
iShankar et"aLll2011[ ). Then the AGN luminosity is 

Xagn = 3.6 X 10^3 erg Mg . (11) 

3.1. AGN spectra 

We adopt the general relativistic, radiatively ef- 
ficient, steady-state th in accretion disk model of 
iNovikov fc Thornd ()1973D . This model describes the flux 
of thermal radiation from t he disk as a function o f radius 
in Equation (5.6.14b) (see lPage &: Thornd 11974 for the 
explicit form of ^) . We assume no radiation from within 
^isco- In addition to the thermal component, AGN 
spectra typically exhibit emission lines, excess infrared 
radiation from dust reprocessing UV emission, and a 
hard X-ray component usually assigned to inverse Comp - 
ton scattering in a hot corona ()Haardt fc Maraschill9"93l) . 
We do not account for these phenomena, but choose our 
observing bands so that their effect is minimal: an ob- 
servation window around 200 eV is low enough so that 
the thermal component dominates over coronal emission, 
but it is higher than helium Lyman absorption and de- 
tector lower energy limits. It is important to note that 
little is known about the geometry of the corona, and 
simultaneously observing a stellar transit in hard X-ray 
might provide feedback to the models. 

We follow t he accretion disk photo sp here model 
described bv iMilosaylievic fc PhinnevI ()2005D and 
iTanaka fc Menoul (|2010D " : we assume that the dominant 
absorption mechanism is the bound-free process, with 
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an opacity that depends on the frequency and tem- 
perature. We assume that the temperature and thus 
the absorption opacity are constant down to an optical 
depth of one (the "bottom") in the photosphere. We 
can calculate the total flux from the absorption and 
electron scattering opacities and the photosphere bottom 
te mperature, using Equat ions (A13-A15) and (A17) 
of iTanaka fc Mmoul pOlOl ). but using the rela t ivistic 
angular frequency given by INovikov fc Thornd (Il973l ) 
instead of Keplerian velocity. However, as the flux 
is known and the temperature is sought for, we have 
to use a simple iterative process to solve this implicit 
equation for the temperature. We find that usually 
electron scattering dominates the total opacity, but 
in the hottest parts of the accretion disk, absorption 
takes over. Note that the functional form of the specific 
flux differs from a blackbody spectrum, because the 
absorption opacity depends on frequency. For simplicity, 
we assume that the emerging radiation is isotropic in 
the frame comoving with the accretion disk. 

Given the specific intensity of the accretion disk as 
a function of radius and frequency, the observed spec- 
trum is determined by Doppler shift, gravitational red- 
shift, and gravitational lensing. To account for these cf- 
fec ts, we apply the tran sfer function method as described 
bylCunninghamI ()1975D and implemented bv lSpeith et al.l 
(|1995f R. We calculate the radiative efficiency as a func- 
tion of spin as described by Shapiro ct al. (1983) to con- 
vert luminosity to mass accretion rate, which is the input 
parameter of this code. We fix the inclination angle 
(the angle between the spin axis and the line of sight) 
at 60°, so that cosi? = 0.5. The calculated spectrum 
of the accretion disk is shown in Figure [1] for a non- 
spinning black hole (with dimensionless spin a = 0) in 
the top panel, and for a highly spinning black hole with 
a prograde disk with a conservative value a = 0.9 in the 
bottom panel. The frequency is in the source rest-frame 
accounting for gravitational redshift, but a possible cos- 
mological redshift for distant AGNs is not considered. 
The specific luminosity value displayed here is what an 
isotropic source would have to have in order to produce 
the same flux as the AGN does at this specific inclination. 

Transit depth is defined as the blocked flux to out-of- 
transit flux ratio in a given band. Therefore the transit 
depth is between zero and one: zero if the transiting ob- 
ject does not cover any part of the accretion disk; one if 
the object completely blocks radiation (in which case it is 
called an occultation or eclipse). The transit depth varies 
with frequency and the location of the transiting object 
in projection. To be able to efficiently calculate tran- 
sit depths at different positions, we implement a linear 
approximation to the radius-gravitational redshift grid 
generated by the above code to determine these values 
for a light ray parametrized by its projected position far 
from the AGN. Then we employ high order numerical 
approximatior0 to integrate over the stellar disk in the 
projection plane. This gives the blocked specific fiux, 
which we then divide by the total specific flux to obtain 
the narrow-band transit depth. 

^ available at http://www.tat.physi k.uiii-tuebmgen.de/| 
[^ipeith / publ / photon_transtcrtct_d blc. tj 

^ code available from http : / /www.hbloborodko.com/pavel/l 
[?page_id =1879 
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Figure 1. Spectrum and transit depth of an accretion disk around 
a lO'' Mq black hole. Solid red curve: equivalent isotropic lumi- 
nosity per logarithmic bins of frequency. Dashed blue curve: max- 
imum possible narrow-band transit depth caused by a star with 
Rt = 11 -R0 as a function of source frequency. Top (bottom) 
panel presents the case of a Schwarzschild (Kerr) BH with spin 
a = (a = 0.9). 

In addition to AGN spectra. Figure [T] also depicts the 
maximum possible narrow-band transit depth caused by 
an early type main sequence star with i?^, = 11 Rq as 
a function of frequency. The maximum depth of a gen- 
eral transit can be smaller if the star does not transit the 
most luminous part of the projected accretion disk. At 
high frequencies, the most luminous region is more com- 
pact, therefore the transit is deeper. The transit depth 
becomes constant at infrared frequencies less than the 
peak frequency of a black body spectrum with tempera- 
ture of the outer edge of the disk (we set 10'^ Rg in the 
simulations), because here the Rayleigh- Jeans approxi- 
mation applies to every part of the disk. Larger stars 
cause deeper transits (see below). 

3.2. Transit maps 

As stated in Section [U AGN transit observations can 
be used to map distant AGNs along the transit chord 
with unprecedented resolution. We now elaborate on 
this point. Figures [2] and [3] show the transit depth as 
a function the projected position of the transiting ob- 
ject in front of the accretion disk around a non-spinning 
BH for an early type main sequence star [R^, = 11 Rq) 
and a giant star {R^, — 110 Rq), respectively. The 
three panels in both flgures show the transit depth maps 
for different observing frequency: optical g band (top 
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Figure 2. The instantaneous narrow-band transit depth as a func- 
tion of projected position of the transiting object at three different 
frequencies optical {top panel), EUV [middle panel), and soft X- 
rays (bottom panel). A transit hghtcurve corresponds to the values 
along the stellar trajectory in the image plane shown. A non- 
spinning 10^ Mq black hole is at the origin, the observation angle 
is 60° relative to the accretion disk, and transit depths are shown 
for a main sequence O-type star with ij* = 11 Rq = 5.18 Rg. The 
top half of the disk image is more severely distorted by gravita- 
tional lensing since it is farther from the observer than the black 
hole. The left side of the disk rotates towards the observer, thus 
appearing more luminous due to beaming, which results in a deeper 
transit. The transit is deepest if the projected position of the star 
crosses the most luminous regions closest to the SMBH. The black 
hole silhouette and the curved spacetime geometry becomes visible 
in the X-ray transit map. 



panel, A = 477 nm, u = 6.3 x 10^"' Hz), extreme ul- 
traviolet (EUV, middle panel, A = 145 nm = 1450 
A, ly = 2.1 X 10^^ Hz), and soft X-ray (bottom panel, 
hv = 200 eV, = 4.8 X 10^^ Hz, A = 6.2 nm). The 
spatial variations in the transit depth maps imply time- 
variations of the observed AGN brightness as a transit- 
ing object moves across the image. A transit lightcurve 
corresponds to the values along the projected stellar tra- 
jectory in Figures [5] and |31 Conversely, observations of 
the transit lightcurve can be used to reveal the geometry 
of the accretion disk along a line in this image. 
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Figure 3. Same as Figure [J] but for a giant star with i?* = 
110 if,0 = 51.8 -Rg. Here, the star is so large that the transit is 
much deeper, and the image of the accretion disk is smoothed out 
(e.g. the left-right asymmetry of beaming is not apparent). 
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Figure 4. Soft X-ray transit depth map for an AGN with a Kerr- 
BH spin a = 0.9 for i?* = 11 Rq, same as the bottom panel of 
Figure[2]but for a spinning SMBH. The details of the accretion disk 
are smoothed since the ISCO radius in this case is smaller than the 
stellar radius. The accretion disk image is further modified because 
of the Kerr geometry. 
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The maximum possible resolution of such a recon- 
structed image is set by the size of the transiting object 
and the spatial variations of the disk brightness. Since 
the disk temperature increases inwards, the emission at 
higher frequencies is confined to the innermost regions, 
implying that transit measurements at higher frequencies 
can give a higher resolution image (see different panels 
in Figures [2] and [3]). Transits of smaller objects also yield 
a higher resolution image. However, transits of smaller 
objects are less deep and hence more difficult to detect. 

The left-right asymmetry in the figure is due to dif- 
ferent Doppler shifts for regions of the disk moving to- 
wards or away from the observer. The spacetime geom- 
etry leaves an imprint on the image, the top half of the 
disk image is distorted by gravitational lensing close to 
the BH. Remarkably, the black hole silhouette (i.e., the 
lack of emission within the ISCO) becomes directly visi- 
ble in the X-ray image of a transiting O star (see Figure [2] 
bottom panel). 

The transit depth map is also sensitive to the space- 
time geometry both directly through gravitational lens- 
ing and indirectly through the change in the ISCO ra- 
dius. Figure |4] shows the soft X-ray transit map of a 
Novikov-Thorne accretion disk around a Kerr BH with 
spin a = 0.9 (c.f. bottom panel of Figure [5]). In this 
case, the ISCO is smaller than the radius of the transit- 
ing object and the BH silhouette does not appear in the 
image. 

3.3. Transit lightcurves 

To get a handle on the characteristic transit duration, 
let us consider the timescale for the center of a star to 
cross a disk of radius i?AGN centered on the most lu- 
minous part of the accretion disk image for a given fre- 
quency. We choose the value of i?AGN based on the tran- 
sit depth map, depending on the desired transit depth. 
Typically i?AGN ~ 5-1000 i?g. 

The transit duration is 



At = 2il'AGN 



= 4 hours 



GAfsMBH 

Ragn 



10 Ra 



1 mpc 



(12) 



The transit duration is of the order of hours for massive 
main sequence stars closest to the AGN, and weeks for 
the closest giants further out in the circumnuclear star 
cluster (see Table [5] below) . 

Figure [S] shows multiband transit lightcurves of an ac- 
cretion disk due to stars of radius i?* = 11 (top and mid- 
dle panel) and 110 Rq (bottom panel). The horizontal 
axis shows the horizontal position of the transiting star in 
the image plane of Figures [IHH Physical distance in the 
projection plane is displayed in Rg units on the bottom 
axis, and the top axis shows time in hours for a star with 
minimum orbital radius as given in Table [TJ Note that 
further stars will cause a similar lightcurve, only scaled in 
time. The dimensionless black hole spin is zero on the top 
and bottom panels and 0.9 in the middle panel. The pro- 
jected stellar trajectory is horizontal with ?/proj = 1.5 R^, 
which approximately passes through the hottest part of 
the accretion disk. The black hole is behind the origin. 
The curves show the transit depth for the same three 
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Figure 5. Lightcurves of a transit: transit depth on logarithmie 
scale as a function of the horizontal position of the transiting star 
in front of the AGN in units of Rg (bottom axis), or time in case 
r = 5 mpc (top axis). Stellar radius is i?* = 11 Rq = 5.2 i?g 
(top and middle panels) and 110 Rq = 52 i?,g (bottom panel). 
The SMBH spin is a = (top and bottom panels) and 0.9 (middle 
panel). The BH silhouette is larger than the star in the top panel 
only. 
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frequencies as the transit depth maps in Figures [2HS1 
Note that these graphs are different from the usual plane- 
tary transit lightcurves showing flux, but plotting transit 
depth is more convenient when using a logarithmic scale. 
A larger value corresponds to a deeper transit, that is, a 
larger photometric signature. Figure [5] is consistent with 
the expectations on the frequency dependence described 
above: at higher frequencies, the AGN is more compact, 
implying shorter and deeper transits. The sharp features 
near the center of the transit on the top and middle pan- 
els are due to the inner edge of the disk at the ISCO, and 
the left-right asymmetry is mostly due to Doppler shift 
of the approaching and receding parts of the disk. How- 
ever, a transiting giant filters out this spatial feature due 
to its size, as seen on the bottom panel. The black hole 
silhouette is clearly visible in the top panel where the 
ISCO radius is larger than the transiting object. Com- 
paring the top two panels, we infer that if the mass of the 
SMBH is known, observing a transit light curve allows 
us to put an upper limit on the spin. 

We conclude that for an accretion disk around a 
10^ Mq SMBH, ~ 0.1% photometric accuracy is required 
in the optical, ^ 1% accuracy in the extreme ultraviolet, 
and ~ 10% accuracy in soft X-ray to detect a transit due 
to a 11 Rq star. In order to detect a transit due to a 
110 Rq giant, ^ 1% accuracy is sufficient in the optical, 
and ^ 10% accuracy in extreme ultraviolet to X-ray (see 
Section 15.11 for a comparison with intrinsic variability) . 

For comparison, we ran simulations for two different 
SMBH masses: 10^ and lO"^ Mq. Table [2] shows the 
maximum possible transit depth in each case for differ- 
ent sizes of transiting stars. The transit depth for a fixed 
stellar size decreases with the increasing spatial extent 
of the accretion disk around black holes with increas- 
ing mass, making the transit feature more prominent for 
smaller BH masses. However, AGNs with less massive 
SMBHs are also intrinsically fainter, and 10^ Mq AGNs 
also have a smaller local number density which makes 
low mass AGN transit observations more challenging (see 
Section [5] for a discussion). 

4. TRANSIT RATES 

In this section, we estimate the expected rates of stellar 
transits in AGNs. Recall that i?AGN denotes the radius 
of the circular area that the center of the star has to 
transit in the projection plane for a given transit depth. 
Let Pinst('') denote the probability of a single star on a 
circular orbit with radius r transiting this circular area at 
a given instance, and let Pever(''') denote the probability 
that the orbit is aligned so that the star transits this 
circular area at some point during its orbit. By geometric 
arguments, these probabilities are 



Pinstir)- 

-Pcvcr(^) ■ 



1 ^AGN 



2i?AGN X 27rr i^AG^ 



47^7-2 



(13) 
(14) 



For a single AGN, the expected value of the number of 
transits at a given instance, iVjnst, the expected value of 
the transit rate (number of transits observed per time for 
asymptotically long observations), N, and the expected 
value of the number of stars on orbits such that they 
ever transit, Never, can be calculated by integrating for 



the entire stellar population: 

iVinst, = / 47rr2n(r)Pi„st(r)dr , 



N = 



47rr^rt(r) 



Pever (^) 



dr , 



Nr.- 



47rr^n(r)Povcr(7-)dr , 



(15) 
(16) 
(17) 



where Torb is the Keplerian orbital period. 

To interpret these probability indicators, we have to 
understand the relationship between the expected value 

N of the transit rate and the expected value iVovcr of 
the number of stars that ever transit a given AGN. If 
we monitor a single target for which A^ovor ^ 1, then 
the actual transit rate is zero with probability 1 — Never 
and N /Never with a small probability Never- (The prob- 
ability of multiple ever transiting stars in a single AGN 
is negligible in this case.) However, when monitoring a 
large number n ^ 1/iVovcr of identical targets, these ef- 
fects average out: the observed total transit rate is « nN 
and there are ~ nNcver stars causing transits among all 
targets in total. 

We substitute Equations into Equations (fT5HT7| . 
and use 1 < a < 3 and rmin ^ ?'max (justified by Table[l]) 
to obtain 



iVinst = 3.8 X 10"" 1000° 



a — 1 



-M. 



3-7(3-a) 



M, 



30 M, 



^ f Ragn \ ^ 

vio rJ 



N = 



1 pc 

1000"-^ 6-^^^ m| 

93 yr - i ^ 



1 



1 mpc 

-7(3-a) 



a - 



/ Ragn \ f rp 
[lORj \lpc 

3 — a 



^30Mo 
N;^>^ = 0.032 X 1000""2 fj 

M^ 



1 mpc 



-M. 



2-7(3-a) 



30 M, 



^ / Ragn \ f 
\IQRJ Vlpc 



6 

Q-3 



X (18) 



(19) 



(20) 



1 mpc 



2-a 



N^=^ = 0.032 X 6(3 - a)Af6"^^^""^ x 



30 



/ Ragn \ f rp 
[lORj \lpc 



In- 



N^<^ = 0.032 X 1000"^^ b- 



-M, 



2-7(3-a) 



a 



(21) 



(22) 



A/. 
30 M© 

ra>2 



/ Ragn\ 
\lORgJ \l pc 



Q-3 



1 mpc 



Here N^v^'^, N"^^, and iV^cr^ denote the value of 
TVover for the corresponding population density expo- 
nents. Stars on close orbits transit more frequently, and 
they dominate TVinst and N . However, Never is dominated 
by stars on close or wide orbits for a > 2 and a < 2, re- 
spectively. This is determined by the exponent of r in 
the integrand of Equations p^HTT)) . 

Let us substitute specific values to make numerical pre- 
dictions. For concreteness, we study the cases of AGNs 
with mass 10^, 10^, and lO'^ Mq, and spin a = and 0.9. 
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Table 2 

Transit probabilities and duration for different bands and transit depths for a single AGN 



M 




model 


X 
(nm) 


maximum 
transit 

dentil^ 


trctnsitr 


n 

^AGN 
(R "I 




1 /N 

(y-i 


AT 

'ever 


At (r ■ ) 

\i mm / 


i-^^ \i max J 

I"; 


10^ 





O star 


477 


0.038 




180 


Q C 
O.O 


X 


iU 


280 


0.009 


89 


750 








145 


0.076 


10-2 


235 


D.o 


X 


in — 6 
iU 


210 


0.012 


120 


980 








6.2 


~ 1 


~ 1 


45 




X 


in-7 

iU 


1 100 


0.002 


22 


190 


10^ 


0.9 


O star 


477 


0.049 


10-2 


200 


4.7 


X 


lU 


250 


0.010 


99 


830 








145 


0.12 


10-2 


240 


6.8 


X 


lU 


210 


0.012 


120 


1000 








6.2 


~ 1 


~ 1 


40 


i.y 


X 


iU 


1300 


0.002 


20 


170 


10^ 





red giant 


477 


0.63 


10-1 


900 


1.2 


X 


10—"* 


1000 


0.11 


830 


20 000 








145 


0.78 


10-1 


750 


O.D 


X 


1 n — 6 
iU 


1200 


0.091 


690 


17 000 








6.2 


0.74 


10-1 


530 




X 


1 n — 6 
iU 


1700 


0.064 


490 


12 000 


10^ 


0.9 


red giant 


477 


0.68 


10-1 


850 


1. i 


X 


in — 5 
iU 


1 100 


0.10 


780 


19 000 








145 


0.88 


10-1 


750 


O.D 


X 


in — 6 
iU 


1200 


0.091 


690 


17 000 








6.2 


~ 1 


10-1 


500 


3.8 


X 


1 n — 6 
10 


1800 


0.060 


460 


11000 


10^ 





O star 


477 


0.0012 


10-^ 


6 


5.6 


X 


lU 


320 


0.010 


1.5 


7.9 








145 


0.0047 


10-3 


50 




X 


in — 5 
iU 


39 


0.083 


13 


66 








6.2 


0.33 


10-1 


9 


1 o 


X 


1 n — 6 
iU 


220 


0.015 


2.3 


12 


10^ 


0.9 


O star 


477 


0.0021 


10-3 


14 


Q n 
O.D 


X 


1 n — 6 
iU 


140 


0.023 


3.5 


18 








145 


0.0089 


10-3 


70 


1.0 


X 


in — 5 
iU 


28 


0.13 


18 


92 








6.2 


0.56 


10-1 


8 


9.9 


X 


in-7 
lU 


240 


0.012 


2.0 


11 


10^ 





red giant 


477 


0.047 


10-2 


220 




X 


in — 5 
iU 


170 


0.78 


98 


2 900 








145 


0.15 


10-1 


55 


O.O 


X 


1 n — 6 
iU 


700 


0.20 


25 


710 








6.2 


~ 1 


~ 1 


40 


Z.O 


X 


1 n — 6 
iU 


960 


0.14 


18 


520 


10^ 


0.9 


red giant 


477 


0.064 


10-2 


220 


8.4 


X 


10-5 


170 


0.78 


98 


2 900 








145 


0.23 


10-1 


75 


9.8 


X 


10-*^ 


510 


0.27 


34 


980 








6.2 


~ 1 


~ 1 


40 


2.8 


X 


10-*^ 


960 


0.14 


18 


520 


10^ 





O star 


477 


0.000027 


10-5 


17 


5.8 


X 


lo-* 


4.4 


0.95 


2.2 


7.1 








145 


0.00016 


10-1 


5 


5.1 


X 


10-5 


15 


0.28 


0.64 


2.1 








6.2 


0.0096 


10-3 


8 


1.3 


X 


10-1 


9 


0.45 


1.0 


3.3 


10^ 


0.9 


O star 


477 


0.000055 


10-5 


35 


2.5 


X 


10-3 


2.2 


2.0 


4.5 


15 








145 


0.00032 


10-* 


12 


2.9 


X 


10-4 


6.3 


0.67 


1.5 


5.0 








6.2 


0.025 


10-3 


6 


7.3 


X 


10-5 


13 


0.34 


0.77 


2.5 


10^ 





red giant 


477 


0.0020 


10-3 


17 


5.8 


X 


10-5 


96 


1.8 


3.7 


130 








145 


0.011 


10-2 


3 


1.8 


X 


10-*^ 


540 


0.32 


0.65 


23 








6.2 


0.55 


10-1 


12 


2.9 


X 


10-5 


140 


1.3 


2.6 


91 




0.9 


red giant 


477 


0.0036 


10-3 


35 


2.4 


X 


10-1 


47 


3.7 


7.6 


270 








145 


0.021 


10-2 


12 


2.9 


X 


10-5 


140 


1.3 


2.6 


91 








6.2 


0.83 


10-1 


8 


1.3 


X 


10-5 


200 


0.84 


1.7 


61 



^ The maximum transit depth is determined by model parameters: SMBH mass and spin, cireumnuclear stellar population, 
and observing wavelength. 

^ This transit depth is chosen to be less than the maximum transit depth. Subsequent columns show values dependent on this 
parameter. The same transit depth is chosen for corresponding a — and a — 0.9 cases to demonstrate how little probability 
estimates depend on spin. 



Table [2] presents these examples. Each line of this table is 
generated by selecting a SMBH mass and spin, a stellar 
population model, and an observing wavelength. We de- 
termine the maximum possible transit depth by running 
our simulation with these input parameters. Then we 
choose a lower value as transit depth threshold, and use 
the calculated transit depth map (similar to Figures [2]-[3l) 
to identify the contour enclosing the area that a star has 
to transit to produce a lightcurve signature of at least 
this depth. We approximate this region by an ellipse, 
and take -Ragn to be the radius of a circular disk of the 
same area. These values can be substituted into Equa- 
tions (IT2|l and p8H22| to calculate the transit duration 
and to obtain probability estimates corresponding to the 
chosen transit depth treshold. 

Table [2] shows that one cannot reasonably expect to 
detect transits when observing only a single or even a 
fistful of AGNs: typical transit rate is one every few 
hundred years (except for the very shallow transits for 
MsMBH = 10^ Mq which may happen more frequently 



than one in 10 years). However, these rate estimates are 
sensitive to poorly constrained parameters such as r^nin, 
b, and the number density exponent a. Indeed, if we 
neglected the limits imposed by stellar collisions and set 
rijiin to the tidal disruption radius, r„iin would decrease 
by a factor of ~ 30 or ~ 4 (see Table [T|) , and the instan- 
taneous transit rate would increase a factor of ~ 200 or 
^ 3, for O stars and giants, respectively. The event rates 
in a flattene d, edge-on oriented star cluster may also be 
much larger (jSeth et al.ll2"008H [2009I ). 

5. PROSPECTS FOR AGN TRANSIT OBSERVATIONS 

In the previous section we established that it is neces- 
sary to monitor a large number of AGNs to confidently 
detect transit events. This can be done either by a tar- 
geted survey, or by monitoring a large region on the sky. 
In this section, we calculate the number density of suit- 
able AGNs in the local universe, then discuss the feasibil- 
ity of observing AGN transits with specific instruments. 
Note that there are other instruments potentially able to 
detect a transit, and archival data of previous observa- 
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Table 3 

Local active black hole density in three decades 





SMBH mass range 




number density'' 








(Mpc-^) 


104.5 


Mq < MsMBH < lO^-S 


Mq 


5 X 10"^ 


105.5 


M@ < MsMBH < 106-5 


Mq 


7 X lO-*^ 


106.5 


M© < MsMBH < lO''-'^ 


Mq 


1 X 10-5 



^ from the lognormal fit given by IGreene fc Hoi 1 1200 71 . 
[2009h 



tions might already contain transits. 

5.1. AGN density and variability 

The transit hghtcurves are sensitive to the assumptions 
on the AGN, e.g. on its mass and Eddington ratio. Larger 
mass means larger luminosity in optical to soft X-ray 
(but the disk is cooler, therefore less luminous in hard 
X-ray). Larger mass also means more extended accretion 
disk, that is, shallower transits. To quantify these effects, 
we consider three decades of magnitude ranges, centered 
on AfsMBH = 10^, 10^, and lO'^ Mq. We integrate the 
lognormal fit of the local act i ve black ho le mass function 
determined by Greene fc Hoi (j2007l[2009D to estimate the 
number density of AGNs in each mass range. Our results 
are given in Table [3] 

For concreteness, we assume that all of these SMBHs 
are highly spinning and have a prograde accretion disk. 
However, Table[2]shows that probabilities depend weakly 
on the spin, therefore this assumption does not influ- 
ence our estimates. (An exception is the case of X-ray 
observations, because accretion disk luminosity in this 
frequency depends strongly on spin, resulting in more 
potential targets, and thus more observable events down 
to a given luminosity limit.) 

We assume the probabilities of the inclination cos i) = 
0.5 case for all AGNs (see Section [5] for a discussion of 
this assumption). 

A transit can only be detected in the lightcurve if 
the AGN variability amplitude on the timescale of in- 
terest is small enough compared to the transit depth. 
IMeusinger et al.l ()2011[ ) investigate a sample of over 9000 
quasars in the SDSS sample between z = 0.2 and 3, and 
find that 93%, 97%, 93%, 87%, and 37% of them are vari- 
able in the u, g, r, i, and z band, respectively. However, 
the variability is dominated by timescales of months to 
years, much larger than the AGN transit timescale of 
hours. A useful measure of the time dependence of the 
AGN variability is the structure function (SF), which 
essentially measures the RMS magnitude difference as a 
function of time lag r betwee n magnitude me a surem ents. 
Based on the SDSS sample, iMacLeod et al.l (|2012f ) find 
that for T < 3 year, SF = 0.02 mag (r/10 day)"-^'*. 
These observations are consistent wit h the damped ran - 
dom walk model of AGN variability (jKellv et al.ll2009D . 
Substituting a timescale of 1 hour, 1 day, and 1 week into 
this relation yields an average AGN variability of 0.0018, 
0.007, and 0.017 magnitudes in the optical g band, which 
corresponds to a variability of 0.17%, 0.7%, and 1.6% 
over these timescales, respectively. Comparing this to 
the AGN transit lightcurve on Figure |5| for transiting gi- 
ants shows that the transit can cause a larger change in 
luminosity than the average intrinsic optical AGN vari- 



ability, provided that the observation cadence is at most 
a few days. Thus we conclude that red giant transits 
may be detectable even in typically variable AGN. How- 
ever, AGN variability may be a limitation for detecting 
transiting O stars in the optical bands, where the transit 
depth is much smaller. 

AGN variability i s more significant on the transit 
timescale in X-ra ys ()Gonzalez-Martin fc VaughanI 120121 : 
iPonti et"aI1l2012D . However, in these bands, the tran- 
sit may be a nearly complete eclipse (Figure jSJ , making 
them detectable regardless of variability. Indeed, tran- 
sits of broad line clouds have alr eady b een detected in 
X-ravs (iMcKernan fc Yaaoob|[T998;.Maiolino et al.l l2010l: 
iRisaliti et al.l l2011bl: iLohfink et al.l I2012D . Broad hne 
clouds are expected to be i nuch more densely distributed 
aroun d AGNs then stars (|Arav et al.l 119971 : ILaor et ahl 
120061 ). However, their transit shapes are different from 
those of stars due to their cometary shape, with high 
column density head s followed by lower column tails 
(|Maiolino et al.l [20101 ). Some of these transiting clouds 
may represent t he irradiated envelopes of circu mnuclear 
"bloated stars" (Edwards' 1980; ' PenstonI IT981) on very 
close orbits (Lohfink ct al. 2012). 

5.2. Ground-based optical instruments 

Pan-STARRS (|Kaiser et al.|[200l IMot ) is an optical 
and NIR survey project, consisting of four units with 
1.8 m primary mirror diameter and a field of view of 
7 deg^ each. The photometric precision is « 1% in most 
bands per ~ 40 s exposure. The first telescope has been 
operating in science mode since 2010. Once all four units 
are online, the system will survey the night sky once 
every week with a 5a detection limit down to r ~ 24. 
We use the conservative magnitude limit g pa 23 for 1% 
photometric accuracy per single exposure, and calculate 
for ten years of operation of all four units. 

An AGN with A/smbh -^10*^ M©, spin a = 0.9, and 
Eddington ratio 0.3 has g = 23 magnitude if observed 
from a luminosity distance of ~ 330 Mpc. (For compar- 
ison, a similar AGN with a = would have the same 
g magnitude from 390 Mpc.) Using cosmocalc.pjQ 
to calculate the comoving volume up to this luminos- 
ity distance, and using the AGN number density in Ta- 
ble 121 we find that there are approximately 600 AGNs 
in the given mass range to this distance observable from 
a single geographic location. Main sequence O stars do 
not cause deep enough transits in the optical bands for 
10^ Mq SMBHs, therefore we focus on transiting red gi- 
ants. Assuming all targets host red giants, there will be 
^ 35 transits in ten years according to the transit rate 
N — 1/(170 yr) given in TablejH When such an event oc- 
curs, the AGN is on the night sky with probability « 1/2. 
The transit lasts longer than a week, therefore there will 
be at least one observation during the transit. Assuming 
0.75 of the dark hours have photometric conditions at 
the site, we estimate that the Pan-STARRS survey can 
detect ~ 10 stellar transits in ~ 10^ Mq AGNs. (If all 
targets had a = 0, the expected number of transits would 
be - 20.) 

A similar calculation for highly spinning 10^ Mq AGNs 
gives an approximate 80 Mpc distance limit for a g = 23 

* by Tom Aldcroft, |http: //cxc. harvard ■ edu/contrib/cosmocalc/| 
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magni t ude. However, in the catalog of iGreene fc Hoi 
()2007l |2009[ ). there is only one active black hole with 
MsMBH < 10^'^ Mq closer than this distance, there- 
fore a Pan-STARRS detection of a transit event in the 
lightcurve of such an object is unlikely in the standard 
all-sky survey mode. Finally, there are approximately 
80 000 AGNs with masses ~ 10^ Mq within 1.65 Gpc 
corresponding to g w 23 mag. We do not expect 1% deep 
transits in this case. The transit rate at a 10"^ photomet- 
ric level is 1/(47 yr), therefore the Pan-STARRS survey 
would observe ^ 6 000 such events during 10 years of op- 
eration. Note however that detecting these transits may 
be prohibitively difficult due to systematics of ground- 
based observations, and also due to AGN variability. 

Future ground-based optical surveys of larger col- 
lecting area might have an even better chance to de- 
tect transits of 10^ Mq active SMBHs. For exam- 
ple, LSST is a survey telescope with an equivalent pri- 
mary mirror diameter of 6.68 m, and a field of view 
of 9.6 deg^, cu rrently in design and development phase 
()LSST Science Collaborations ct al. 2009) . It is planned 
to have an observing strategy similar to that of Pan- 
STARRS, but with 3.4 times the collecting area of the 
four Pan-STARRS units in total, it can survey an ap- 
proximately 3.4'^/^ w 6 times larger volume, increasing 
the transit detection expectation accordingly. Therefore 
we estimate that LSST may detect ~ 100 transits per 
decade for Msmbh ~ 10*^ Mq. Also, since LSST is 
planned to have more frequent visits than Pan-STARRS, 
the intrinsic variability will have a smaller amplitude be- 
tween subsequent observations. 

Since many survey telescopes devote a certain fraction 
of their time to deeper surveys of smaller areas, it is 
worth estimating how this changes the expected number 
of transit detections. Recall that the four Pan-STARRS 
units will be able to scan the entire night sky approxi- 
mately once a week in five filters with « 40 second expo- 
sures. Now assume a single Pan-STARRS class instru- 
ment spends all available time on surveying an n times 
smaller area in g filter only, with m visits per week to 
fight AGN variability, thus exposing n/m times longer. 
As long as the observations are photon noise limited and 
cosmological effects are negligible, this strategy increases 
the surveyed volume, and also the expected number of 
transit detections, by a factor of n^/^rn~^/^. For exam- 
ple, daily visits (m = 7) of an n w 1300 smaller area 
(a few fields) would mean two hour exposures per point- 
ing, and would double the expected transit detections 
around Msmbh ^ 10^ Mq AGNs to - Z-i per year, 
with increased robustness against intrinsic variability. 

5.3. Kepler observations 

The Kepler satellite carries out almost continuous pho- 
tometric measurements in the bandpass of 420 to 900 nm. 
It achieves « 10"'' photometric a ccuracy per 30 mi nute 
e xposure on Kp = 1 3 dwa rf stars ()Koch et al.ll2010[ ). 

IMushotzkv et al.' (2011") repo rt ob servations of four 
AGNs, and Edelson & MalkinI ()2012[ ) identify 13 more 
in the Kepler field. However, these quasars are at red- 
shifts between z = 0.028 and 0.625, therefore it is not 
likely that either of them has low enough SMBH mass 
that transits due to O stars or red giants could be ob- 
served by Kepler. 



5.4. Space X-ray instruments 

The Chandra X-ray observatory has been operating 
in orbit since 1999. Its relevant instrument is the High 
Resolution Camera (HRC) has a 31 arcmin by 31 ar- 
cmin field of view, se nsitive from 0.1 keV to 10 keV 
(|WeisskoDf et al.l [200^. 

The X-ray Multi-Mirror Mission (XMM) Newton satel- 
lite is another proposal instrument, operating since 2000. 
It is equipped with th ree imaging cameras sensitive from 
0.2 k eV to 12 keV HTurner et all l2Q0ll: iStriider et al] 
120011) , with approximately 700 arcmin^ field of view each. 

The Spektr-RG satellite is scheduled to launch in 2013 
at earliest. Its eROSITA X-ray telescope system has 
approximately twice the effective area o f a single in- 
strument on XMM-Newton below 2 keV (jPredehl et al.l 
|2010[) . This observatory is planned to carry out a sur- 
vey consisting of eight full scans of the sky in four years, 
with a mean total exposure time of 2 ks for each region 
(Brunner ct al. 2011). 

The Extreme Physics Explorer (jGarcia "eTalllMll ) IS a 
mission concept designed specifically to study accretion 
disks around SMBHs. It is proposed to have more than 
an order of magnitude larger effective area than current 
missions, thus capable of observing targets farther away. 

As an example, we estimate the expected number of 
transits detected during a single 100 ks observing cam- 
paign with the HRC-I instrument on Chandra. We con- 
sider observations in the energy range 0.1 keV to 0.4 keV, 
logarithmically centered on the energy 0.2 keV for which 
we calculated transit probabilities. Based on Section [3TT1 
the integrated luminosity in this range is 9 X 10^^ erg/s, 
8 X 10''^ erg/s, and 6 x lO''^ erg/s, and the photon index is 
-0.2, 0.4, and 1.6, for AGNs with a = 0.9 and mass 10^ 
10^, and 10'' Mq, respectively. (Accretion disks around 
non-spinning black holes are up to a factor of four less 
luminous in this energy range.) 

Let us divide the observation time into ten bins of 10 ks 
each for temporal resolution. If we want the relative pho- 
ton noise to be 0.1 in each bin, we need a count rate of at 
least 0.01 photons per second. According to the PIMMS 
count rate calculatoi0, this is the case up to a luminosity 
distance of w 290, 890, and 2 300 Mpc for SMBH masses 
of 10^, 10^, and 10^ Mq, respectively. This corresponds 
to a redshift z w 0.06, 0.2, and 0.4, which is small enough 
to neglect in terms of flux change for a simple estimate. 
There are « 40, « 12 000, and « 180 000 AGNs in the 
corresponding mass ranges, out to this luminosity dis- 
tance, out of which « 0.0003, w 0.08, and « 1 fall in the 
Chandra field of view on average. Therefore we can only 
observe a single or at most a few targets in one pointing. 
Since the probability of detecting a transit when observ- 
ing a single target for this much time is negligible based 
on the transit rate value of Table [H such a short cam- 
paign is not suitable for stellar AGN transit discovery. 

Now let us consider the planned eROSITA survey. Sup- 
pose the 2 ks exposure time is evenly distributed among 
eight visits of 250 s on each object. We have the same 
0.1 relative photon noise per visit out to approximately 
200 Mpc, 700 Mpc, and 2 Gpc, for AGNs with 10^ 10^, 
and 10^ Mq, respectively. On the entire sky, this means 
approximately 20, 6 000, and 100 000 AGNs. For such 

® |http : //cxc ■ harvard ■ edu/toolkit/plmms • j sp | 
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short observations, the expected number of transits sam- 
pled is the product of the number of AGNs, the number 
of visits per AGN, and the instantaneous probabihty of 
a transit, iVinst- Based on the vahies in Table EJ if we 
assume O stars around each target, we expect 3 x 10~^ 
almost total echpses for AGNs with mass ~ 10^ M©, 
and 0.05 transits of depth 0.1 for AGNs with 10^ Mq. 
In case of SMBH mass around lO'' Mq, transits due to 
O stars do not reach the depth of 0.1, therefore they are 
not hkely to be detected in the inherent variability. If 
we assume a b = 0.01 mass fraction of giant stars in a 
typical AGN, we expect approximately ~ 10~^, 0.1, and 
10 transits of depth ~ 1, ~ 1, and 0.1 detected for SMBH 
masses ^ 10^, 10^, and 10'' Mq, respectively. 

We conclude that discovering stellar transits in X-ray 
with proposal instruments is unlikely, whereas it may be 
possible with survey instruments like eROSITA for AGNs 
with lO^-lO"^ Mq. 

6. DISCUSSION 

In this section we highlight the uncertain points in our 
argument, not only to be able to properly interpret our 
predictions, but also to understand the implications of 
future AGN transit detections or non-detections. 

We based our stellar models in galactic nuclei on ob- 
servations of the center of the Galaxy in [J2l assuming 
that this is a good representation of circumnuclear star 
clusters in active galaxies, and on theoretical dynamical 
models. We set up simplified stellar population models 
to estimate transit rates. Our estimates of transit prob- 
abilities depend strongly on the assumed value of the 
stellar distribution exponent a: the 1000""^ terms in 
Equations (fT8 l) -(f22 l) are due to the three order of magni- 
tude difference in the distance scales of r^in and tq. As 
a consequence, steeper stellar distributions feature larger 
transit probabilities, consistent with the statement that 
most probabilistic measures are dominated by closeby 
stars for moderate values of a. 

Evaluated at a = 3, Equation (|3]) and thus Equations 
(fT8l - [22|) yield zero. This artifact originates in the as- 
sumption that the number of stars within rmin is negli- 
gible, which is not true for this value of a. To study the 
case of a w 3, we recalculate N without this assumption: 



N = 



0.093 yr a - i ^ 



M, 
30 Mf 



(23) 



1 mpc 



/ Ragn\ 

[lo rJ 



Figure [S] displays this formula for as a function of a. 
We plot two cases: transit rate with a g-band threshold 
depth of lO^'^ for O stars (i?AGN = 14), and with a g- 
band threshold depth of 10~^ for red giants or bloated 
stars with = 110 Rq, M^ = 1.5 Mq (i?AGN = 220). 
We fix Me = 1, a = 0.9, and b = 0.01 for both cases. The 
plots show us that TV is a steeply increasing function of 
a up to and over the value 3, as expected. Comparing 
the plot to the values given in Table [2] at a — 2.5 in the 
first case {N — 1/140 yr) and at a = 1.75 in the second 
case {N = 1/170 yr) tells us that our approximation in 
Equation ((3)) for these values of a is valid. 
Our method to determine the number of stars within 
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Figure 6. Transit rate N for a single target AGN as a function 
of stellar density exponent a, with Mg = 1, a = 0.9, and b = 0.01. 
Solid red curve corresponds to O stars at a g-band transit depth 
threshold of 10^'^, dashed green curve corresponds to red giants or 
bloated stars at a g-band transit depth threshold of 10~^. 

the sphere of influence based on the M-a relation is over- 
simplified. The actual number of stars depends on for- 
mation, relaxation, and disruption rate and history. Fur- 
thermore, the presence of stellar remnants may dilute 
the stellar population within the sphere of influence for 
a fixed total mass, decreasing the mass fraction b - and 
therefore the number - of stars that would potentially 
cause detectable transits. 

At large orbital radii, we expect microlensing to domi- 
nate over the light obscuration due to the transiting star. 
We did not investigate this phenomenon in detail but 
conservatively ignored all transits due to stars on wide 
orbits where microlensing may become significant. As 
the transit rate is dominated by closeby stars, microlens- 
ing is not a limiting factor of transit observability after 
all. Also, we did not account for apsidal and Lense- 
Thirring precession, which may move stars in and out of 
transiting orbits. 

AGN transits due to broad line clouds are expected 
to be much more freque nt than the stellar t r ansits 
calculated here (se e e.g. IMcKernan fc YaaoobI 119981: 
Bianchi et al. 2009: jMaiolino et all l2010t iRisaliti et al l 
2007fc^2009atbi,,201ibl!a!). Transits of clouds may be some- 



what different than s tellar transits due to their cometary 
shape (|Maiolino et al.ii20 10'). and typically have a shorter 
duration due to their proximity to the SMBH and larger 
velocity (typically At ^ 1 hr Mq if orbiting at tblr ^ 

1000 Rg). 

Note that we have assumed that the transiting stars 
move on orbits much wider than the accretion disk. This 
assumption may be violated for stars on very close or- 
bits (?■ < 10^ Rg = Mq 4.8 mpc), where stars cross- 
ing the disk may get captured by hydrodynamic drag 
(|Karas fc Subdl200l . 

We mentioned the possibility of supermassive stars 
possibly forming in AGN accretion disks, but do not 
have information on their occurrence rates. Furthermore, 
Wolf-Rayet stars exhibit strong stellar winds, which 
might form a n optically thick region of radius ^100 Rq, 
(Figure 5 in iCrowtheil 120071 and references therein) ; 
OH/IR stars can form dus t clouds much larger than 
the e ntire accretion disk fe.g. lBednnlll987l : lKemper et al.l 
l2002i) : and bloated stars with large irradiated envelopes 
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(|Edwardslll980t IPenstonlll988l) may be present near the 
AGN. These objects can potentially cause much deeper 
transits or even eclipses. 

In Section 13.11 we investigated the projected ther- 
mal radiation structure of an accretion disk with spe- 
cific parameters, assuming radiatively efficient accretion, 
modelling radiative transfer in the accretion disk photo- 
sphere, and accounting for light propagation in the Kerr 
metric of the SMBH. We set the Eddington ratio to 0.25 , 
a value based on the analyses of iKoUmeier et aP ()2006f ) 
and[Shankar ct al. (20 1 ^. T h is val ue is consistent with 
the findings of iGreene fc Hd (|2007t ). However, keep in 
mind that magnitude limited samples are biased towards 
larger Eddington ratios. Also note that the Eddington 
ratio can in fact vary b y an order of m agnitude in either 
direction (e.g. Kauffma nn fc Heckma n 2009). A larger 
accretion rate would increase the total luminosity and 
the peak frequency with little effect on the size of the 
accretion disk. Therefore such a disk would be much 
brighter in X-ray than one with lower Eddington ratio, 
while they would exhibit transits of similar depth. 

We also assumed prograde disk alignment for spinning 
SMBHs, but noted that the transit probabilities do not 
depend strongly on spin. We fixed the inclination at 
cosi? = 0.5, which is the mean value for an isotropic 
distribution, therefore the predicted transit probabilities 
are typical of all inclinations. However, note that larger 
inclination results in a thinner image of the accretion 
disk, therefore deeper transits. A disk with thickness 
H / R « 0.05 observed from a nearly edge-on orientation 
results in a '-^ 10 times deeper transit, leading to a 10'^ 
times larger detectable volume for photon noise limited 
surveys. On the other hand, there are selection effects at 
both inclination extrema: a coplanar torus might obscure 
thermal emission from edge-on AGNs, whereas jets along 
the rotational axis might contaminate the lightcurve of 
face-on AGNs and prevent transits from being detected. 
These are likely to confine the inclination distribution 
of AGNs with observable thermal radiation closer to the 
average values we use in our model. 

Intrinsic AGN variability will pose a challenge to iden- 
tifying AGN transits. Fortunately, optical AGN vari- 
ability is small on the timescales of days, and thus does 
not rule out the possibility of transit observations for gi- 
ant stars. With Pan-STARRS observations, however, the 
variability level at typical observation cadence is compa- 
rable to the transit depth. Transits might not be detected 
if the weekly observations miss the deepest part. 

Simultaneous multiband observation campaigns can 
help distinguish variability from transit signiture, as the 
transit lightcurve is predicted to have a different shape 
in different frequency bands. Also, transit lightcurves 
could be contaminated by the flux reflected by clouds 
surrounding the AGN. Multi-wavelength campaigns fol- 
lowing AGN for these types of events may give informa- 
tion about the reflecting fraction, which could allow to 
constrain the covering fraction of material, the location 
and relative sizes of the reflecting regions. Future de- 
tectors could distinguish the reflected component using 
polarimetry. 

7. CONCLUSIONS 

In this paper, we presented simple estimates to study 
the prospects for detecting stellar AGN transits. We have 



shown that such observations would offer a novel possi- 
bility to image the accretion disk in distant AGN with 
unprecedented accuracy. For example, the black hole sil- 
houette (i.e., the lack of emission within the ISCO) may 
be resolved using the lightcurve of an AGN transit due to 
an 0-type main sequence star. These observations probe 
the accretion disk and the space-time geometry around 
black holes, and in particular, they are sensitive to the 
black hole spin. AGN stellar transit event rates offer 
information about the circumnuclear stellar cluster. 

We predict that the Pan-STARRS survey could detect 
10 ~ 20 stellar transits, and LSST may detect ^ 100 by 
repeated photometric observations of ~ 10^ Mq AGNs. 
We estimate that stellar transit detections in X-rays are 
not likely with individual campaigns on proposal in- 
struments, but we expect a few possible detections in 
~ 10^ Mq AGNs with X-ray surveys like eROSITA. 
Note that these rate estimates do not include transits 
by Compton thick clouds, which are observed to be com- 
mon in X-rays. The transit rate corresponding to clouds 
could also be much more frequent in optical surveys. 

However, these probability estimates are very sensitive 
to parameters which are based on theoretical arguments. 
One of these parameters is the inner radius cutoff of the 
nuclear stellar cluster, rmin, which is set by stellar col- 
lisions. We have shown that if stellar collisions did not 
deplete the innermost regions, the less stringent limita- 
tion due to tidal disruption would imply ~ 200 times 
larger transit rates for O stars. Such scenarios may 
be possible if the effective stellar size increases during 
close approach to the AGN, le ading to large i r radiated 
envel opes ( "bloated stars" , see lEdwardsl 119801 : iPenstonI 
|1988() . AGN transit observations could constrain these 
parameters and reflne circumnuclear stellar population 
models in distant galaxies. 
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